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We find all the exact eigenstates and eigenvalues of a spin-1/2 model on square lattice: H — 
^^o'^Zi ^i+x^i+x+y^i+y- show that the ground states for g < and g > Q have different 
quantum orders described by Z2A and Z2B projective symmetry groups. The phase transition at 
g = Q represents a new kind of phase transitions that changes quantum orders but not symmetry. 
Both the Z2A and Z2B states are described by Z2 lattice gauge theories at low energies. They have 
robust topologically degenerate ground states and gapless edge excitations. 
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Introduction 

We used to believe that all phases of matter are de- 
scribed by Landau's symmetry breaking theory. ||l|, || The 
symmetry and the related order parameters dominate our 
understanding of phases and phases transitions for over 
50 years. In this respect, the fractional quantum Hall 
(FQH) states discovered in 1982 j|, || opened a new chap- 
ter in condensed matter physics. The theory of phases 
and phase transitions entered into a new era. This is be- 
cause all different FQH states have the same symmetry 
and hence cannot be described by the Landau's theory. 
In 1989, it was realized that FQH states, having a ro- 
bust topological degeneracy, contain a completely new 
kind of order - topological order. Q A whole new theory 
was developed to describe the topological orders in FQH 
liquids. (For a review, see Ref. |^.) 

The Landau's theory was developed for classical statis- 
tical systems which are described by positive probability 
distribution functions of infinite variables. FQH states 
are described by their ground state wave functions which 
are complex functions of infinite variables. Thus it is not 
surprising that FQH states contain addition structures 
(or a new kind of orders) that cannot be described by 
broken symmetries and the Landau's theory. From this 
point of view, we see that any quantum states may con- 
tain new kind of orders that are beyond symmetry char- 
acterization. Such kind of orders was studied in Ref. 
and was called quantum order. Since we cannot use or- 
der parameter to describe quantum orders, a new mathe- 
matical object - projective symmetry group (PSG) - was 
introduced to characterize them. The topological or- 
der is a special case of quantum order - a quantum order 
with a finite energy gap. 

One may ask why do we need to introduce a new con- 
cept quantum order? What use can it have? To answer 
such a question, we would like to ask why do we need 
the concept of symmetry breaking? Is the symmetry 
breaking description of classical order useful? Symmetry 
breaking is useful because (A) it leads to a classification 
of crystal orders, and (B) it determines the structure of 
low energy excitations without the needs to know the de- 
tails of a system. P, p[] The quantum order and its PSG 



description are useful in same sense: (A) PSG can classify 
different quantum states that have the same symmetry, 0] 
and (B) quantum orders determine the structure of low 
energy excitations without the needs to know the details 
of a system. 0, |ll| The main difference between classi- 
cal orders and quantum orders is that classical orders pro- 
duce and protect gapless Nambu-Goldstone modes[^, ||] 
which is a bosonic excitation, while quantum orders can 
produce and protect gapless gauge bosons and gapless 
fermions. Fermion excitations (with gauge charge) can 
even appear in pure bosonic models as long as the boson 
ground state has a proper quantum order. ^ ^ 

Those amazing properties of quantum orders could 
fundamentally change our views on the universe and its 
elementary building blocks. The believed "elementary" 
particles, such as photons, electrons, etc , may not be 
elementary after all. Our vacuum may be a bosonic 
state with a non-trivial quantum order where the "el- 
ementary" gauge bosons and the "elementary" fermions 
actually appear as the collective excitations above the 
quantum ordered ground state. It may be the quantum 
order that protects the lightness of those "elementary" 
particles whose masses are 10^° below the natural mass 
- the Plank mass. Those conjectures are not just wild 
guesses. A quantum ordered state for a lattice spin model 
has been constructed! 10 which reproduces a complete 



*URL: littp: //dao .mit . edu/'wen 



QED with light, electrons, protons, atoms, ... 

The concept of topological/quantum order is also use- 
ful in the field of quantum computation. People has been 
designing different kinds of quantum entangled states to 
perform different computing tasks. When number of 
qubits becomes larger and larger, it is more and more 
difficult to understand the pattern of quantum entangle- 
ments. One needs a theory to characterize different quan- 
tum entanglements in many-qubit systems. The theory 
of topological/quantum order[g[ Q is just such a theory. 
In fact topological/quantum orders can be viewed as pat- 
terns of quantum entanglements and gauge bosons the 
fluctuations of quantum entanglements. Also the robust 
topological degeneracy in topological ordered states dis- 
covered in ReLJ a] can be used in fault-tolerant quantum 
computation. llS f 

It is hard to convince people about the usefulness of 
quantum orders when their very existence is in doubt. 
However, a growing list of soluble or quasi soluble 
models |jl|, |l|, ^ ||, |l^, Q indicates that topolog- 
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ical/quantum order do exist beyond FQH states. In 
particular, Kitaev has constructed exactly soluble spin 
models that realize both topological orders {ie quan- 
tum orders with finite energy gap) and gapless quantum 
orders. 

In this paper, we study an exact soluble spin- 1/2 model 
on square lattice: H = 16gj:, Sy Sf^^S^^^^^Sf^^. We 
find that the ground states for g < and g > have 
the same symmetry but different quantum orders. The 
PSG's for those quantum ordered states are identified. 
The phase transition at g = represents a new kind of 
phase transitions that changes quantum orders but not 
symmetry. We show that the projective construction that 
is used to construct quantum ordered ground states] 12 
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1| , |22| | not only gives us exact ground states for 
our model, but also all the exact excited states. Through 
this soluble model, we hope to put quantum order and 
its PSG description on a firm ground. 

We would like to mention that the above spin- 1/2 
model, having one spin per unit cell, is different from Ki- 
taev's exact soluble spin-1/2 models on the links of square 
lattice and on the sites of honeycomb lattice (which have 
two spins per unit cell).|jl^, Q However, the g < 
version of the above model corresponds to the low en- 
ergy sector of Kitaev's honeycomb lattice model in the 
Jz » Jx,Jy limit. |3) 

Quantum orders in spin-1/2 and hard-core-boson 
models 

In this section, we are going to give a brief and general 
review of the PSG description of quantum order. Read- 
ers who are interested in the exact soluble model can go 
directly to the next section. Let us consider a spin-1/2 
system on a square lattice. Such a system can be viewed 
as a hard-core-boson model if we identify | |) state as 
zero-boson state |0) and | t) state as one-boson state |1). 
In the follow we will use the boson picture to describe 
our model. 

To construct quantum ordered (or entangled) many- 
boson wave functions, we will use projective con- 
struction. We first introduce a "mean-field" fermion 



Hamiltonian:[^ 



mean — ^ ^ 
(ij) 



4,iXl^ i^J-J + ^\,d^^\j + h.c) (1) 



PSG, we first discuss two types of transformations. The 
first type is SU (2) gauge transformation 

{^i^Xio^mj) ^ {G{i)^i,G{i)xijG\j),G{i)inijG^{j)) 

(4) 

where G{i) S SU{2). We note that the physical bo- 
son wave function <I>*^^^^''''^') ({i„}) is invariant under the 
above SU{2) gauge transformations. The second type is 
the usual symmetry transformation, such as the transla- 
tions Tx'. i ^ i — X, Ty: i i — ij. A generic trans- 
formation is a combination of the above two types, say 
GT^iXij) = G(i)xi-i,j-xGt(i). The PSG for an ansatz 
iXijiVij) is formed by all the transformations that leave 
the ansatz invariant. 

Every PSG contains a special subgroup, which is called 
the invariant gauge group (IGG). An IGG is formed 
by pure gauge transformations that leave the ansatz 
unchanged IGG = {G\ Xij = G{i)xijG^ {j),Tjij 
G{i)T]ijG'^{j)}. One can show that PSG, IGG, and the 
symmetry group (SG) of the many-boson wave function 
are related: PSG /IGG = SG.^^ 

Different quantum orders in the ground states of our 
boson system are characterized by different PSG's. In 
the following we will concentrate on the simplest kind 
of quantum orders whose PSG has a 100=^2. We will 
call those quantum states Z2 quantum states. We would 
like to ask how many different Z2 quantum states are 
there that have translation symmetry. According to our 
PSG characterization of quantum orders, the above phys- 
ical question becomes the following mathematical ques- 
tion: how many different PSG's are there that satisfy 
PSG/Z2 ^translation symmetry group. This problem 
has been solved in Ref. . The answer is 2 for 2D square 
lattice. Both PSG's are generated by three elements 
{GxTx,GyTy,Gg\, whcrc Gg is a pure gauge transfor- 
mation that generates the Z2 IGG: IGG = {1, Gg}. The 
gauge transformations in the three generators for the first 
Z2 PSG are given by 
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(5) 



Such a PSG wiU be called a Z2A PSG. The quantum 
states characterized by Z2A PSG will be called Z2A 
quantum states. For the second Z2 PSG, we have 



where I,J = 1,2. We will use Xij a-nd Vij to denote 
the 2x2 complex matrices whose elements are xlj and 



rjl-. Let l^-^^ean 



be the ground state of the above 
free fermion Hamiltonian, then a many-body boson wave 
function can be obtained 



'J>{il,i2...) = 



in 



6(i„)|*^e^j;I'-^) (2) 



where 



b{i) = Vi,iV'2 



(3) 



According to Ref. , the quantum order in the boson 
wave function $'^^*3'''*J'^({i„}) can be (partially) charac- 
terized by projective symmetry group (PSG). To define 



Ggii) = -1, G,(i) = l, Gy{i)^i-iy 



(6) 



Such a PSG will be called a Z2B PSG. 

If wc increase the symmetry of the boson wave func- 
tion, there can be more different quantum orders. A 
classification of quantum orders for spin-1/2 system on a 
square lattice is given in Ref. [@] where hundreds of dif- 
ferent quantum orders with the translation, parity, and 
time-reversal symmetries were found. 

In the following we will study an exact soluble model 
whose ground states realize some of the constructed 
quantum orders. Our model has the following nice prop- 
erty: the projective construction Eq. (||) gives us all the 
energy eigenstates for some proper choices of Xij and rjij . 
All the energy eigenvalues can be calculated exactly. 
Exact soluble models on 2D square lattice 
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Our construction is motivated by Kitaev's construction 
of soluble spin-1/2 models on honeycomb lattice. The 
key step in both constructions is to find a system of com- 
muting operators. Let U^j = ^JU^■\j, where i,j label 
lattice sites, a is an integer index, Ufj is an n x n matrix 
satisfying {U^-Y = -U^^, and Xj = (Ai,i, A2,i, A„,i) 
is a n-component Majorana fermion operator satisfying 
{Xa,i,Xbj} — 25ab5ij- We require that all U^^ to com- 
mute with each other: \tjf ■ ■ ] — 0, which can be 
satisfied iff 

[/?.[/?>. =0 ■ U- ■ = (If^ ■ ■ V 

[/? . . =0 (71 

where ii, i2, is are all different. 

Let {|1), |4)} be a basis of a 4 dimensional real lin- 
ear space. Then the following Uij on a square lattice 
UiA+s. = |1)(3|, C/i,i_i = -|3)(1|, Ui^i+^ = |2)(4|, and 
Ui.i-y = -|4)(2| form a solution of Eq. (0). We find that 

Ui^i+x — Ai^jAs^i+x, Ui^i+y = \2A^iA+y (8) 

form a commuting set of operators. 

After obtaining a commuting set of operators, we can 
easily see that the following Hamiltonian 

H .9 ^ ^ 1 — Ui_i-^tji-^_i^tji^_i^tji^_i (9) 

i 

commutes with all the Uij^s, where ii — i + i2 — 
i + X + y, and is = i + y. We will call Fi a Z2 flux 
operator. Let \sij) be the common eigenstate of Uij with 
eigenvalue Sij. Since (Uij)'^ = —1, Sij satisfies sij = ±i 
and Sij — —Sji. \sij) is also an energy eigenstate with 
energy 

F — ^ ^ Fi , Fi — Si^i^Sij^^i2Si2^i^Si^^i (-^0) 
i 

Let us discuss the Hilbert space within which the above 
H acts. On each site, we group Ai^2,3.4 into two fermion 
operators 

V'l^i = Ai,i iAa^i, ^2,i = A2,i -1- iA4,i (11) 

"01,2 generate a four dimensional Hilbert space on each 
site. Let us assume the 2D square lattice to have Ng 
lattice sites and a periodic boundary condition in both 
directions. Since there are total of 2^^= different choices 
of Sij (two choices for each link), the states \sij) exhaust 
all the 4^= states in the Hilbert space. Thus the com- 
mon eigenstates of Uij is not degenerate and the above 
approach allows us to obtain all the eigenstates and eigen- 
values of the H. 

We note that the Hamiltonian H can only change the 
fermion number on each site by an even number. Thus 
the H acts within a subspace which has an even num- 
ber of fermions on each site. The subspace has only two 
states per site. When defined on the subspace, H ac- 
tually describes a spin-1/2 or a hard-core boson system 
under the operator mapping Eq. (H). The subspace is 



formed by states that are invariant under local Z2 gauge 
transformations: ipn G{i)ipii, G{i) = ±1. We will 
call those states physical states and call the subspace the 
physical Hilbert space. 

Since Uij do not act within the physical Hilbert space, 
they do not have definite values for physical states. How- 
ever, we note that the Z2 flux operator Fi act within 
the physical Hilbert space. The Z2 flux operators com- 
mute with each other and the if is a function of Fi. 
To obtain the common eigenstates of the Z2 flux opera- 
tors in the physical Hilbert space, we note that the Z2 
flux operators are invariant under the Z2 gauge trans- 
formation: ipii — > G{i)ipii. A Z2 gauge transforma- 
tion changes one eigenstate \sij) to another eigenstate 

(Hi G'f'')|sij} oc \G{i)sijG{j)), where Ni is the fermion 
number operator at site i. We will call those two eigen- 
states gauge equivalent. The common eigenstates of the 
Z2 flux operators within the physical Hilbert space can 
be obtained by summing over all the gauge equivalent 
eigenstates with equal amplitude: 

G i 

(Note \sij)phy can also be viewed as the projection of 
\sij) onto the physical Hilbert space.) 

Let us count the physical states that can be con- 
structed this way, again assuming a periodic boundary 
condition in both directions. We note that the terms in 
the above sum can be grouped into pairs, where the Z2 
gauge transformations in a pair differs by a uniform Z2 
gauge transformation G{i) = — 1. We see that the sum 
of the two terms in a pair is zero if the total number of 
the fermions Nf =J2-^i odd. Thus only states \sij) 
with even number of fermions leads to physical states 
in the above construction. For states with even num- 
ber of fermions, there are only 2^" /2 distinct terms in 
the above sum. Thus each physical eigenstate comes 
from 2^' /2 gauge equivalent eigenstates. Since there are 
A^' /2 states with even number of fermions, we find the 
number of physical states is 2''^= , which is the dimension 
of the physical Hilbert space. Thus we can obtain all 
the eigenstates and eigenvalues of the H in the physical 
Hilbert space from our construction. 

Let us introduce a notion of Z2 flux configuration. 
Among all the possible Sij's, we can use the Z2 gauge 
transformation Sij G(i)sijG{j) to define an equiva- 
lence relation. A Z2 flux configuration is then an equiv- 
alence class under the Z2 gauge transformation. From 
the above picture, we see that for every choice of Z2 
flux configuration Fi = ±1, we can choose a Sij that 
reproduces the Z2 flux on each plaquette. If the state 
\sij) has even numbers of fermions, it will leads to a 
physical eigenstate (see Eq. (p^)). The energy of the 
physical eigenstate is given by Eq. (|o|). The explicit 
many-boson wave function of the eigenstate is given by 

*J'({*«}) = (0|nn^(^«)l*mean) wherC |*mean) is the 

ground state of 

Hmean = ^ {sijUij + h.c)j (13) 

(ii) 
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Here we see that all the eigenstates of our model can be 
obtained from the projective construction Eq. (^). 
Physical properties 

In terms of the hard-core-boson operator Eq. (||), the 
Hamiltonian Eq. (|) of our soluble model has a form H = 
~9Y.i{b-b'^)i{b+V^)i+x[b-h'^)i+x+y{h+h^)i+y- In terms 
of spin-1/2 operator = 6 + 6^ and — i{b — b^), the 
Hamiltonian has a form 



H 



9 



bx+y 



(14) 



In the bulk, the above model corresponds to the following 
simple Ising model Hi sing = 9j2i'''i- However, the two 
models are different for finite systems and for systems 
with edges (see discussions below). 

When g < 0, the ground state of our model is given 
by Z2 flux configuration Fi — 1. To produce such a flux, 
we can choose Sj^i+i — Si i^y ~ i. In this case, Eq. dl^ ) 
becomes Eq. (|l|) with -rii^i+^ = Xi,i+x = 1 + r"' and 
—Vi,i+y = Xi,i+y — 1 — T^. The PSG for the above 
ansatz turns out to be the Z2A PSG in Eq. (||). Thus 
the ground state for g < is a Z2A state. 

When g > 0, the ground state is given by configura- 
tion Fi = —1 which can be produced by (— )*''Si.i+x = 
Si,i+y ^ i- The ansatz now has a form — ?7i.i+x = 
Xi,i+A = (-)*"(! + t"") and -rji^i+j. = Xi.i+y = I - t\ 
Its PSG is the Z2B PSG in Eq. (|). Thus the ground 
state for 5 > is a Z2B state. 

Both the Z2A and the Z2B states has translation T^.y, 
parity Pxy ■ {ix,iy) {iy,ix)i and time reversal (since 
Xij and "qij are real) symmetries. The low energy excita- 
tions in both states are Z2 vortices generated by flipping 
the signs of an even numbers of Fis. Those Z2 vortex 
excitations have a finite energy gap A = 2\g\. The low 



energy sector of our model is identical to the low energy 
sector of a Z2 lattice gauge theory. However, our model 
is not equivalent to a Z2 lattice gauge theory. This is 
because a Z2 lattice gauge theory has 4 x 2^=/2 states 
on a torus while our model has 2^= states. 

Due to the low energy Z2 gauge structure, both the 
Z2A and the Z2B states have four degenerate ground 
states on an even by even lattice with periodic boundary 
condition. The degeneracy is topological and is protected 
by the topological order in the two states. The degen- 
eracy is robust against arbitrary perturbations. This is 
because perturbations are Z2 gauge invariant physical 
operators which cannot break the low energy Z2 gauge 
structure that produces the degeneracy. ||, |l2| 

On an even by even lattice, the ground state energy per 
site is given by —\g\. The singularity at g = implies a 
phase transition between two states with the same sym- 
metry and the same ground state degeneracy! Thus the 
transition is a new type of continuous transitions that 
only changes the quantum orders (from Z2A to Z2B). 
The transition is continuous since gap vanishes at (7 = 0. 
Although the Z2A and the Z2B states share many com- 
mon properties on an even by even lattice, the two states 
are quite different on an odd by odd lattice. On an odd 
by odd lattice, the Z2A state has an energy — IgjA^s and 
a 2 fold degeneracy, while the Z2B state, containing a 
single Z2 vortex to satisfy the constraint Fj = 1, has 
an energy — I^KA^s — 2) and a 2Ns fold degeneracy. On 
lattice with edges in (1,0) and/or (0,1) directions, our 
model has '-^ 2^=''9=/^ gapless edge states, where N^dge is 
the number of edge sites. 
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